Trapping and cooling fermionic atoms into the Mott and Neel states 
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We perform a theoretical study of a fermionic gas with two hyperfme states confined to an optical 
lattice. We derive a generic state diagram as a function of interaction strength, particle number, 
and confining potential. We discuss the central density, the double occupancy and their derivatives 
as probes for the Mott state, connecting our findings to the recent experiment of Jordens et al. [lj. 
Using entropic arguments we compare two different strategies to reach the antiferromagnetic state 
in the presence of a trapping potential. 
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Remarkable experimental progress has been made in 
handling and controlling quantum degenerate atomic 
gases. The confinement of these gases within optical lat- 
tices, and the use of Feshbach resonances allow for the in- 
vestigation of strongly correlated quantum phases [2|] . An 
outstanding example is the experimental study of the su- 
pcrfluid to Mott insulator transition in bosonic gases [3j . 
A Mott insulating state is characterized by the suppres- 
sion of density fluctuations. Very recently, the observa- 
tion of an incompressible Mott state has been reported 
in a fermionic gas using two hyperfme states of 40 K [1(. 
A future major step will be the observation of the anti- 
ferromagnetic ordering of the two hyperfme components 
in the lattice. However, whereas an incompressible Mott 
state is reached when the temperature is low enough as 
compared to the Mott gap, stabilizing a magnetically or- 
dered state requires cooling down to a lower temperature 
scale set by the superexchange inter-site coupling. 

In this article, we address theoretically several issues 
which have immediate relevance to those experimental 
efforts. We establish a state diagram of cold fermions 
subject to an optical lattice and a trapping potential, as 
a function of the number of atoms and coupling strength. 
We discuss how to detect different states, by considering 
observables such as the fraction of doubly occupied sites, 
the local occupancy and the compressibility at the cen- 
ter of the trap. Finally, we outline two possible strategies 
for adiabatically cooling the system into a long-range or- 
dered antiferromagnet. 

A quantum degenerate mixture of two hyperfme states 
of fermionic atoms can be described by a Hubbard-type 
Hamiltonian [4|, |5(: 
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Here, a =t, | labels the two hyperfme states, j the sites of 
a three-dimensional (3D) cubic lattice, and (j, j') neigh- 
boring sites. Tj is the distance between j and the center, 
d the lattice spacing, and pq the chemical potential. We 



consider a mixture with equal number N/2 of f and J, 
atoms. Experimentally the interaction parameter U and 
the hopping parameter J can be tuned independently 
changing the optical lattice height and using a Fesh- 
bach resonance. Note that the description by a one-band 
model is only valid for large enough lattice depth, small 
enough scattering length, and low temperatures, i.e. sit- 
uations in which the higher bands and more complicated 
interaction or hopping terms can be neglected [J, [fj[ . 

One important effect of the confining potential is to 
allow for the spatial coexistence of different quantum 
states [7|, |8|. For a ID tube of fermionic atoms, a Mott- 
insulating phase next to a liquid phase was pointed out 
[3, [lO| . More recently, the stability of antiferromagneti- 
cally ordered regions has been investigated in a 2D sys- 
tem [ll|] . Helmes et al. [l2j studied the coexistence of 
the liquid and insulating phases, in 3D, calculating in 
particular the density and momentum profiles and the 
local spectral function. In Fig. [TJ we display a state di- 
agram which summarizes the different possible shapes of 
the density profiles as a function of coupling strength 
and atom number at the temperature fcsT/(6J) =0.1. 
All calculations in this paper were obtained using a local 
density approximation (LDA), in combination with dy- 
namical mean field theory (DMFT) [13j, [21 1 to simulate 
the homogeneous system in the paramagnetic state. 

In order to compare different trapping energies and 
particle numbers, we use the scaled number of atoms 
p = N{V t /6J) 3 / 2 (see Ref. 19] for ID). Within LDA 
and the continuum limit, p entirely determines the shape 
of the density profile at a fixed interaction strength and 
temperature, allowing for a straightforward comparison 
to experimental setups. The state diagram in Fig. [1] dis- 
plays four distinct regimes, labeled by (L), (Mc), (Ms) 
and (B). (L): For small interaction or small density, 
the whole system is in a (Fermi-) liquid phase. There, 
the density profile depends strongly on the trapping 
potential. (Mc): Above a critical interaction strength 
u = U/(6J) > 2 and a characteristic density p > 5 a 
Mott central plateau with n « 1 forms. A liquid phase 
surrounds this plateau. (Ms): When the particle number 
is further increased the pressure exerted by the trapping 
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FIG. 1: (Color online) Left panel: State diagram for fermionic 
atoms as a function of interaction strength u = U/(6J) and 
scaled particle number p = N(Vt/6J) 3 ' 2 . Symbols show re- 
sults of DMFT calculations. The boundaries of the region 
(Mc) with a Mott-insulator in the center of the trap are de- 
fined from the constraint no — 1| < 0.005 on the central 
density. Lines are an analytical continuum approximation at 
T = 0. Right panel: typical density profiles for regions (L), 
(Mc) and (Ms), with the radius R set by n(R) = 0.1. 



potential overwhelms the incompressibility of the Mott 
state. This leads to a state with a central occupancy 
larger than one particle per site surrounded by a Mott 
insulating shell. (B): For even larger particle numbers 
a band-insulating plateau with n = 2 forms, surrounded 
by liquid and possibly Mott-insulating shells. 

To detect the occuring quantum phases we investigated 
mainly two observables: the occupancy at the center of 
the trap no and the fraction of atoms residing on doubly 
occupied sites D = -^ Si( n j'T n ii)- D was measured by 
molecule formation in the experiments of Ref. [l]. The 
experimental observation of the density in a small central 
region is more involved. However, a lot of effort is cur- 
rently devoted to the development of probes with good 
spatial resolution. To make contact with current experi- 
ments, we assume the slow experimental loading process 
[1| into the optical lattice to be adiabatic, and perform 
calculations at a fixed value of the entropy per particle, 
s [22j. The latter is related to the initial temperature 
of the non-interacting fermion gas in the absence of the 
lattice by T Z /T F = s/n 2 + 0[(T,/T F ) 3 }. 

The first panel of Fig. [2] displays n and D as a func- 
tion of particle number at a fixed initial temperature 
TijTp = 1.0/-7T 2 . As we detail later on, a plateau in n vs. 
p provides a good indicator of the formation of an incom- 
pressible Mott insulating region in the center of the trap. 
The three sets in Fig. [2] (first panel) correspond to differ- 
ent regimes of coupling. For weak interactions (u < 1.5), 
no plateau is present: no smoothly increases and satu- 
rates at no ~ 2 corresponding to the regimes (L) and 
(B). At intermediate coupling u > 2.5, a Mott plateau 
starts forming {p ss 5), corresponding to regime (Mc). no 
eventually deviates from unity as p is further increased 
into regime (Ms). With increasing coupling (cf. u > 5), 
the central plateau broadens. From Fig. [3] (first panel), it 
is also clear that there is a qualitative correlation between 




FIG. 2: (Color online) First panel: Occupancy no at the cen- 
ter of the trap and fraction D of atoms residing on doubly 
occupied sites versus the scaled particle number p at con- 
stant initial temperature Ti/TF = s/-k 2 = 1.0/-7T 2 . The cen- 
tral occupancy no and the "compressibilities" dno/dp and 
dD/dp, for u = 5, at two different initial temperatures 
Ti/TF = 1.0/7T 2 , 1.6/7T 2 are shown in the second and the third 
panel, respectively. 



the existence of a Mott central plateau and the behavior 
of D. At weak coupling, D just increases smoothly as p 
is increased. At stronger coupling, D remains very small 
for the low values of p, and starts raising approximately 
at the transition from the central incompressible regime 
(Mc) into regime (Ms). 

We now make this correlation between D and no more 
quantitative ploting in the lower panels of Fig. [2] the 
derivatives dno/dp and dD/dp, for u = 5, at two dif- 
ferent initial temperatures Ti/Tp = 1.0/ir 2 , 1.6/ir 2 . In 
the LDA approximation, dno/dp is directly proportional 
to the local compressibility at the center of the trap 
Ko = dno/dpo, through: pdno/dp = Kq/ '(k). In this ex- 
pression, (k) is an average of the compressibility over the 
whole trap. The average (k) is a non-singular quantity 
since it always has sizeable contributions from the liquid 
wings of the density profile. Hence, dno/dp is a direct 
measure of the local incompressibility and its vanishing 
signals the Mott state. Indeed, it is seen from Fig. [5] 
(second panel) that, for the lowest temperature studied 
Ti/Tp = 1.0/vr 2 , this derivative is very small (< 1CT 4 ) 
inside the (Mc) regime, and increases dramatically as one 
departs from a Mott insulator in the trap center. For the 
higher temperature Ti/Tp = 1.6/-7T 2 (Fig. [21 third panel), 
dno/dp is only of order 10~ 3 in the rather narrow (Mc) 
regime, and increases less steeply as p is increased. In 
both cases, it is apparent from Fig. [2] that, remarkably, 
the dependence of dD/dp on p is very similar to that 
of the 'compressibility' dno/dp, for values of p in the 
(Mc) and (Ms) regimes. The reason for this observation 
is that, despite being a global quantity over the whole 
trap, the double occupancy fraction remains very small 



in the boundary regions due to the low density. Therefore 
it is dominated by the central regions with occupancies 
ri > 1. Hence dD/dp is small in regime (Mc) and in- 
creases sharply when the incompressibility at trap center 
is broken. At smaller values of p, these quantities are 
obviously very different: D is very small due to the low 
filling whereas dn^/dp is finite due to the liquid character 
of the state. Thus, we conclude that D and dD/dp are 
rather good indicators of when the Mott state ceases to 
exist at trap center. This validates the use of these quan- 
tities as a diagnostic for a Mott insulator region, as done 
in Ref. [l(. A closely related observation was made by 
Scarola et al. [14| |. who further supported it by showing 
that the compressibility dn/dp is essentially identical to 
d(D / N) / d p in the homogeneous system for n > 1. 

Comparing directly to the experimental results shown 
in Fig. 2 of Ref. [lj , we find good agreement between the 
upturn of D(p) in the experimental data for U/6J = 4.8 
approximately at p — 20 and our theoretical results for 
C//6J = 5 (cf. Fig. [5]). However, in this regime of 
coupling, the estimated lowest temperature in the ex- 
periments (Ti/Tp ~ 0.15), while being at the onset of 
the incompressible Mott regime, still implies a sizeable 
contribution of thermal excitations within the Mott gap 
(Fig- El third panel). In contrast, at lower temperature 
(second panel) a much flatter plateau and lower value 
of the compressibility are obtained. Obviously, a Mott 
plateau forms at even higher temperature for stronger 
interactions, but then equilibration times may become 
an issue. 

A crucial point in experiments is to estimate the tem- 
perature of the interacting fermionic system in the pres- 
ence of an optical lattice. For a confined non-interacting 
gas the double occupancy fraction has been used exper- 
imentally [151 ] as a thermometer since D(T) decreases 



as the temperature is increased [16J, |17[. We find that 
the decrease of D with temperature persists throughout 
the weak and intermediate interaction regimes (u < 1.5, 
Fig. [3l lower panel) . As in the non- interacting gas, the 
transfer of atoms from the trap center (cf. Fig. up- 
per panel) towards the trap boundaries, due to thermal 
fluctuations, causes this effect. It leads to a decrease 
of D which overwhelms the increase found in the high- 
T or large-C/ regime for the homogeneous system. Also, 
the Pomeranchuk effect [6J, [l8| at low-T and intermediate 
U reinforces this trend. In contrast, for stronger inter- 
actions, D hardly changes with temperature (u < 2.5, 
Fig. [3| and ceases to be useful for thermometry because 
the density profile becomes almost frozen as the incom- 
pressible Mott state is approached. 

One major step after observing the Mott regime will be 
the realization of the antiferromagnetically (AF) ordered 
phase, which occurs at a still lower temperature TV (Neel 
temperature). We now discuss two possible strategies to 
reach this phase, (a) A first possibility is to work at 
intermediate coupling (u ~ 2.5), for which TV is the 
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FIG. 3: (Color online) Upper panel: Central occupancy no 
versus temperature for different interaction strengths. Lower 
panel: double occupancy D versus temperature, (p « 13.3) 



largest. (For possible ways of increasing TV further, see 
[l9(). The drawback of this strategy is that a sizeable 
Mott insulating region can only be stabilized within a 
narrow range of particle numbers (cf. Fig. [I]), (b) The 
second possibility is to work at stronger coupling (e.g. 
u ~ 5) so that a Mott region can be stabilized over a 
much wider range of particle numbers (Fig. [TJ. Since 
Tn oc J 2 /U at large U, this requires cooling the gas 
down to very low temperatures. However, assuming an 
adiabatic evolution, the entropy per particle at which 
the ordering takes place, s(TV), is the key quantity to 
focus on, rather than the absolute temperature TV. In 
Ref. [6j, it was pointed out that, for the homogeneous 
half-filled system, s/j(TV) reaches a finite value in the 
strong-coupling limit, even though TV is very small in 
this limit. Here, we compare the initial temperatures 
Ti/Tf needed to reach the AF ordered state for the two 
strategies (a) and (b), in the presence of the trapping 
potential [23l | . 

Once an incompressible Mott region is formed in the 
center of the trap, the entropy is the sum of two contri- 
butions: one from the liquid wings of the density profile, 
and one from the central Mott region. The entropy of the 
Fermi-liquid fraction is linearly increasing with temper- 
ature, below some characteristic coherence temperature 
(which is larger at low density and smaller for densities 
close to one particle per site). By contrast, the entropy 
in the Mott state remains almost constant upon cooling. 

At intermediate coupling (a), the contribution of the 
liquid wings dominates over that of the central Mott re- 
gion in the whole temperature range above TV (Fig. [H 
left panel). This is because the fraction of atoms in the 
central region is rather small (~ 30%) and TV in this 
regime is high enough such that part of the wings have 
not yet reached Fermi-liquid coherence. Hence, strategy 
(a) can benefit from the high entropic contributions of 
the liquid wings. The lowest temperature displayed in 




FIG. 4: (Color online) Entropy per particle s of the total 
system (filled symbols) and the Mott region (open symbols) 
versus characteristic density p. The lowest T is chosen slightly 
above the Neel temperature and the arrows delimit the den- 
sity range where a central Mott plateau is formed. 



coupling route (a) is somewhat more favorable from the 
cooling point of view. By contrast the stronger coupling 
route (b) is more favorable in terms of the number of 
atoms and extension of the Mott plateau, while the cor- 
responding ordering temperatures may still be reachable. 

Note added. — After the completion of this work also 
Schneider et al. (arXiv:0809.1464) reported the observa- 
tion of an incompressible Mott state. 

We acknowledge fruitful discussions with E. Demler, 
M. Kohl, W. Ketterle, C. Salomon, M. Zwierlein, the 
group of T. Esslinger and I. Bloch. Support was provided 
by the 'Triangle de la Physique', the DARPA-OLE pro- 
gram, ICAM and the Agence Nationale de la Recherche 
under contracts GASCOR and FABIOLA. 



Fig. 0] is close to T N at u = 2.5. The entropy per par- 
ticle which must be reached is of order s ~ 0.67 (taking 
p ~ 5 — 10 where a Mott central region is stabilized). This 
corresponds to an initial temperature Ti/Tp ~ 0.07. 

For stronger interaction (b), the density profile is dom- 
inated by a wide Mott-insulating region, for a large range 
of particle numbers. The entropy of this region gives a 
major contribution to the total entropy already at rela- 
tively large temperatures (~ 50% at fcsT/6J = 0.1 as 
shown in Fig. IU right panel) and it does not change 
much upon cooling. In contrast, the contribution of the 
liquid wings decreases strongly and eventually becomes 
very small (30% at k B T/6J = 0.025). In this regime, 
DMFT yields a value of the entropy per particle of order 
0.6 near Tjy, comparable to strategy (a) above. However, 
DMFT overestimates the entropy of the Mott state. We 
propose a lower bound on the entropy at T/y by neglect- 
ing the contribution of the liquid wings, and considering 
only that of the Mott region. This leads to the estimate 
s = s h (T N )N c /N. Using a typical value N c /N = 2/3 
for the fraction of atoms in the central Mott region, and 
the entropy per particle of the homogeneous Mott state 
at Tjy, s/j(Tjv) ~ In 2/2 (as estimated from fluctuation 
effects in the Heisenberg model [6|, see also [20|), we ob- 
tain s ~ 0.2. This corresponds to an initial cooling down 
to Ti/Tp — 0.02 to reach AF ordering. We view this 
estimate as a pessimistic lower bound. 

To summarize we established a state diagram of cold 
fermions in an optical lattice and a trapping potential. 
We further showed how to detect the Mott-insulating 
state using the central density, the fraction of doubly oc- 
cupied sites and their derivatives. The latter were shown 
to be related to the compressibility at the center of the 
trap. We found that the double occupancy is well suited 
for thermometry in the lattice at weak and intermedi- 
ate coupling, but not at strong coupling. Finally, wc 
compared two strategies for adiabatically cooling into a 
long-range ordered antiferromagnet. The intermediate 
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